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ABSTRACT 



We present extensive tables of the atomic partition function [APF) for nickel ions, Ni I 
- Ni X. Partition functions are given over wide range of temperature, 10^ K < T < 10^ K, 
and lowering of ionization energy (0.001 eV < LIE < 5.0 eV), both taken as independent 
variables. Our APF take into account all energy levels predicted by quantum mechanics, 
including autoionization levels. The tables can be applied for the computations of model 
stellar atmospheres and theoretical spectra over very wide range of spectral classes, from 
the coolest K-M dwarfs up to the hottest main sequence, giant, and white dwarf stars. 
This include also model spectra of supersoft X-ray sources and accretion discs in interacting 
binaries. 

Our tables are available at [tittp : //www . astrouw . edu . pi/ '^acta/ acta . htmlj (Acta 



Astronomica Archive) , and [http : / /draco . uni . opole . pl/Halenka . html 
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1. Introduction 

The knowledge of realiable atomic partition functions (APF) is of extreme importance 
for the determination of ionization balance in astrophysical plasmas. In particular, 
partition functions are coefficients of the Saha-Eggert equations describing ionization states 
of elements in stellar atmospheres (equation of state for plasma in Local Thermodynamic 
Equilibrium). Therefore the knowledge of numerical values of the APF is essential 
for computations of model stellar atmospheres and theoretical spectra, and for correct 
interpretation of intensities of absorption spectral lines and abundance determination of 
elements in stellar atmospheres. 

Given value of the partition function U for a particular ion depends on gas temperature 
T, and the local electron concentration A^e- Temperature enters partition function by its 
definition (e.g. Griem 1964, Drawin and Felenbok 1965, Traving et al. 1966): 

U{T, iV,) = exp{-E,/kT) , (1) 

where the sum is taken over discrete energy levels of statistical weight i and excitation 
energies Ei. In vacuum (A^e = 0) the number of bound levels imax is infinite, and therefore 
the above series always diverges. 

In real plasma, however, interaction between the atom of interest and surrounding free 
electrons and ions (plasma effects) cause, that bound levels of very high excitation energies 
move to continuum and no longer contribute to the partitiion function. Therefore the series 
in Eq. (1) reduces to finite number of terms, and value of U is also finite and is strongly 
dependent on the electron concentration A''^. In general, the larger is the lower is both 
number of bound energy levels i„iax in Eq- (1), and the value of U. 

There exist a large number of papers, which present tables of rather approximate (or 
even schematic) partition functions (APF) for elements, including nickel ions (Drawin and 
Felenbok 1965; Traving et al. 1966; Irwin 1981, for example). The latter paper present 
fitting formulae for APF of Ni I - Ni III, for temperatures T < 16000 K. The widely used 
computer code Tlusty 195 for computations of NLTE model stellar atmospheres (Hubeny 
and Lanz 1992, 1995) contains FORTRAN subroutine computing APF of Ni IV - Ni 
IX by direct summation over all observed energy levels of these ions. However, all these 
partition functions for Ni ions depend only on temperature T, and no level dissolution with 
increasing density is included here. The latter implies, that the set of energy levels was 
not complete there. 

Quality of given tables of partition functions depend on (i) accuracy and completeness 
of energy levels included in Eq. (1), and (ii) realiability of the assumed theory of emitter- 
plasma interactions. Emitter-plasma interactions cause that the series in Eq. (1) is finite. 
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Unfortunately, none of the currently existing theories describe correctly effects of charged 
particles in plasma on the atomic partition functions, cf. also Hummer and Mihalas (1988). 

Taking this into account we have decided to compute and present tables of the partition 
functions taking into account all energy levels predicted by quantum mechanics, including 
also levels lying above the so called normal ionization energy (autoionization levels). 
Moreover, our computations arc more physically correct, since the APF depend on both 
temperature T and electron concentration A^^, through tabulated values of lowering of the 
ionization energy (LIE). The method of APF calculations is briefly described in the 
following Section. Detailed description of this method was presented in a series of papers 
by Halenka and Grabowski (1977, 1984, 1986), Halenka (1988, 1989), and Madej et al. 
(1999). 



2. Computation of the atomic partition functions 

We define the atomic partition function of r-th ionization state by the equation 

Pmax Kp)max Pmax 

cf. Halenka and Grabowski (1977). Here the set {pi) with the indices p and i (ordering 
of levels from the ground level upwards in energy scale) describes an eigenstatc of the 
atom in the r-th ionization state. Index i represents three quantum numbers {nlj) of the 
optical electron, and p represents the quantum state of the atomic core. Index i{p)rnax is 
the number of all bound energy levels, Qp^' and Epl^ denote statistical weight and excitation 
energy of the i-th state, in the sequence based on the p-th parent level. Numbers i{p)max 
result from the inequality 

^^S' < - AS^"^) , (3) 

where Ai?'^'"-' denotes the lowering of the ionization energy LIE, and i?^^ is the ionization 
energy in the p-th level sequence. The latter quantity is equal to the sum 

poo loo p \ / 

The quantity E^''^^^ denotes the energy of the atomic core after ionization, r r + 1. 
Index Pmax is the number of different parent levels which can be realized in given physical 
conditions. The number Pmax results from the inequality similar to Eq. (3), written for 
the (r + l)-th ionization state. Since for a fixed value of p the number k is assigned 
unambigously, then the upper limit of E"^^ for the /c-fold excitation {k — 1,2, . . . , Z — r, 
where Z is the atomic number) can be written as follows 

k 

£;M < y . (5) 

poo — / ^ loo \ / 

s=l 
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Following Eq. (3) we have computed extensive tables of the APF for nickel in the 
10 lowest ionization states, Ni I - Ni X. Excitation energies and statistical weights of the 
"observed" levels available for nickel ions were taken from Kurucz (1994). We have added 
to our partition functions contribution from many energy levels predicted by quantum 
mechanics, but missed in his catalogue. This include also many autoionizing levels. Exact 
description of the method of adding of missing levels is given in Halenka and Grabowski 
(1977, 1984). 

In practical model atmosphere calculations one has to estimate the lowering of 
ionization energy for nickel as function of temperature T and electron concentration 
A^e- Such a relation has to result from the model describing plasma-emitter interaction. 
Unfortunately, none of the existing models is satisfactory enough (Hummer and Mihalas 
1988). We suggest use of the approximate relation between the LIE and N^, 

AX = Ze'/D = 3 X 10-' ZN'J^ T-''^ [eV] (6) 

(Eq. 9-106 of Mihalas, 1978), where D = 4.8 {T/N^Y'^ [cm] is the Debye length in hydrogen 
dominated plasma. However, other relations of this type are also given by Drawin and 
Felenbok (1965). 

3. Results 

As the ilustration of our recommended results. Fig. 1 presents run of APF for Ni V and 
temperatures corresponding to atmospheres of hot white dwarf stars, for various values of 
LIE taken as free parameter (solid lines). Numerical values of APF for Ni V are listed in 
Table 1. For a comparison, we have computed also partition functions taking into account 
so called observed levels (single dashed line). One can easily note, that our recommended 
partition functions are significantly larger than the latter functions. 

Partition functions for nickel are arranged in 10 ASCII tables, where each table 
corresponds to a single ion. Complete set of these data is available from Acta Astronomica 
Archive , [http : //www . astrouw . edu . pl/^acta/acta . htinl| , 
or [http : //draco .uni . opole .pl/Halenka.html . 

Entries of a table are decimal logarithms of the APF. They are tabulated at 62 discrete 
temperatures, spaced at nonequidistant intervals, 10^ <T< lO*" K, and at 9 arbitrarily 
assumed values of the lowering of ionization energy {LIE = 0.001, 0.003, 0.010, 0.030, 
0.100, 0.300, 1.000, 3.000, and 5.000 eV). Both T and LIE points remain identical in aU 
10 tables of APF, to ensure homogeneity of the data. 

We are aware, that some values of T and LIE in our tables do not correspond to 
conditions met in astrophysical plasma. However, the extend of all tables ensures, that 
they cover practically all T and LIE expected in stellar atmospheres of any type, excluding 
atmospheres of known neutron stars (temperatures T > 10^ K). 
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Tables presented in this paper are based on the most complete set of energy levels 
actually available. Moreover, our partition functions are sensitive for plasma interactions, 
i.e. they strongly depend on the lowering of ionization energy, which is expected in plasma. 
Our APF tend to diverge for the lowering of ionization energy approaching zero, which is 
the fundamental property of the partition functions in general. 




Fig. 1 : Run of partition functions of Ni V as function of gas 
temperature T and various parameters LIE. Solid lines represent our 
recommended results, whereas dashed line represents partition function 
computed from the observed levels only. 



Acknowledgements 

We are grateful to R.L. Kurucz for making his CD-ROM No. 22 available for us. JH 
and JM acknowledge support by grant No. 2 P03D 013 19 from the Polish Committee for 
Scientific Research. 



6 



REFERENCES 

Drawin, H.W., and Felenbok, P. 1965, Data for Plasma in Local Thermodynamic 
Equilibrium (Gauthier-Villars, Paris) 

Griem, H.R. 1964, Plasma Spectroscopy (McGraw-Hill Book Co., New York) 

Halenka, J. 1988, Astron. Astrophys. SuppL, 75, 47 

Halenka, J. 1989, Astron. Astrophys. SuppL, 81, 303 

Halenka, J., and Grabowski, B. 1977, Astron. Astrophys., 54, 757 

Halenka, J., and Grabowski, B. 1984, Astron. Astrophys. SuppL, 57, 43 

Halenka, J., and Grabowski, B. 1986, Astron. Astrophys. SuppL, 64, 495 

Hubeny, L, and Lanz, T. 1992, Astron. Astrophys., 262, 501 

Hubeny, L, and Lanz, T. 1995, Astrophys. J, 439, 875 

Hummer, D.G., and Mihalas, D. 1988, Astrophys. J., 331, 794 

Irwin, A.W. 1981, Astrophys. J. SuppL, 45, 621 

Kurucz, R.L. 1994, CD-ROM No. 22 

Madej, J., Halenka, J., and Grabowski, B. 1999, Astron. Astrophys., 343, 531 
Mihalas, D. 1978, Stellar Atmospheres (Freeman, San Francisco) 

Traving, G., Baschek, B., and Holweger, H. 1966, Abh. Hamburger Sternw., VIII, No. 1 



Tabic 1 

Atomic partition functions for Ni V (decimal logarithms) 



Lowering of ionization energy (eV) 
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,6247 


29, 


,5071 


23 


.7701 


19, 


.2717 


17 


.1810 


950000 


56, 


.0295 


51 


.0176 


45, 


.5151 


40, 


,4511 


35, 


.0050 


29, 


,8875 


24 


.1484 


19, 


.6277 


17 


.5211 


1000000 


56, 


.3723 


51 


.3605 


45, 


.8580 


40, 


.7941 


35, 


.3478 


30, 


.2303 


24, 


.4895 


19, 


.9493 


17 


.8290 



